Quantum entanglement and interference from classical statistics 
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Quantum mechanics for a four-state-system is derived from classical statistics. Entanglement, 
interference, the difference between identical fermions or bosons and the unitary time evolution 
find an interpretation within a classical statistical ensemble. Quantum systems are subsystems 
of larger classical statistical systems, which include the environment or the vacuum. They are 
characterized by incomplete statistics in the sense that the classical correlation function cannot be 
used for sequences of measurements in the subsystem. 



A classical description of quantum mechanics has been 
tried many times, partly motivated by conceptual difRcul- 
ties of understanding simple quantum phenomena, as the 
Einstein- Rosen-Podolski (EPR) paradox In particu- 
lar, it was hoped that a deterministic description in terms 
of "hidden variables" may be possible. In present days, 
however, it is common belief that the phenomenon of en- 
tanglement, which is the center-piece of the EPR-paradox, 
is genuinely of a quantum mechanical nature and differen- 
tiates quantum statistics from classical statistics. Entan- 
glement is the key stone for quantum computing and 
the basis of theoretical considerations of the foundations 
of quantum mechanics |3j]- Beautiful experiments employ 
entanglement for teleportation 

One may take the attitude that the fundamental descrip- 
tion of the world must be of a probabilistic nature [5] ; with 
deterministic laws arising only in limiting (albeit quite gen- 
uine) cases. In this context we propose here that quantum 
statistics can be obtained from classical statistics. On one 
side it is not a deterministic "hidden variable theory" , on 
the other side it does not need any concept beyond a clas- 
sical statistical ensemble. The differences between classi- 
cal statistics and quantum statistics are only apparent - 
they are due to particularities of an "incomplete" statisti- 
cal description of subsystems 0]. No additional concepts 
are needed beyond the classical statistical ensemble with a 
probability distribution for classical states, fixed values of 
classical observables in these states, and appropriate con- 
ditional probabilities for the outcome of sequences of mea- 
surements. 

Several attempts have shown quantum mechanical fea- 
tures in classical statistical ensembles [1] , : but a full im- 
plementation of all features of quantum mechanics within 
classical statistics has not been achieved previously. In a 
sense, the functional integral representation of quantum 
field theory is already quite close to classical statistics Q. 
The euclidean version after analytic continuation is often 
given by a standard classical statistical ensemble with a 
positive probability distribution. Unfortunately, this nec- 
essary positivity of the probability is lost in the formulation 
in Minkowski space with real time - a complex phase factor 
replaces the euclidean probability distribution. 

If a classical statistical description of quantum mechan- 
ics is possible, an explicit construction should be possible 
for the simplest systems with two or four states. Since four 
states can already describe an entangled state between two 



two-state subsystems, we concentrate in this note on this 
case. We explicitly present a classical statistical ensem- 
ble which realizes all laws of quantum mechanics for the 
four-state system, and demonstrate how entanglement and 
interference arise in this classical statistical setting. 

A classical statistical ensemble which describes four-state 
quantum mechanics typically involves infinitely many clas- 
sical states r, with probabilities Pr > 0, J^tPt ~ ^- Such 
an ensemble describes the quantum system together with 
its environment or the vacuum. The quantum system can 
be regarded as an isolated subsystem, where the mean- 
ing of "isolation" will be discussed later. Only a small 
part of the information contained in the probability distri- 
bution {pt} is needed for a description of the subsystem, 
while the remaining part characterizes the environment or 
the vacuum. More precisely, a pure state of the subsys- 
tem can be described by a set of real numbers fk, which 
correspond to expectation values of observables that can 
be measured in the subsystem, without invoking detailed 
knowledge of the environment. We assume that all prop- 
erties of the subsystem can be described by the set {fk}- 
For four-state quantum systems fifteen numbers fk (which 
obey additional constraints) will be needed. Obviously, 
the specification of fifteen expectation values is only a very 
small part of the information contained in the infinity of 
probabilities Pr- While all fk can be computed from {pr}, 
the inverse is not possible since many different probability 
distributions {pr} lead to the same set {fk}- 

Probabilistic observables 

We concentrate on possible measurements that can only 
resolve two bits. (In a quantum language this corresponds 
to two spins that can only have the values up or down.) 
For any individual measurement, the measurement-device 
or apparatus can only take the values +1 or —1 for bit 1 and 
the same for bit 2. In total there are four possible outcomes 
of an individual measurement, i.e. (++),(H — ),( — h) and 

( ). We describe such measurements by probabilistic 

observables A that are characterized by the probabilities 
w[^\fk) and w[^\fk) = 1 — w^\fk) to find a value -1-1 
or —1 in any given state of the subsystem. These prob- 
abilities only depend on fk- In contrast, for any classical 
state r, the observable has a fixed value Ar = ±1, such 
that = 1 for all states. The probabilities Wj^'' to find 
values A — ±1 are related to the expectation value of A in 
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a given state of the subsystem {fk} 

{A){h) = Y^PrAr = wi^\fk) - W^^\fk) 



(1) 



We note that = w^_^>{fk) + w':^\fk) - 1, as it 

should be for an observable that can only take the val- 
ues ±1 for any individual measurement. (See ref. fS] for 
more details about probabilistic observables and [10] for 
the option that probabilistic observables can be used as 
the fundamental statistical definition of observables.) We 
assume that ^ is a "system observable" in the sense that 
{A) depends only on and not on properties of the envi- 
ronment. 

We focus first on three such "two-level observables", 
namely Ti for the measurement of bit 1, T2 for the mea- 
surement of bit 2, and T3 for the measurement of the prod- 
uct of bit 1 and bit 2. Denoting by w++,w^ ,w |- and 

w the probabilities to measure in the ensemble the out- 
comes (+-I-), (H — ), ( — h) and ( ), one has 



(Ti) = W++ + - W-+ - w — 

{T2) = W++ - w+- + w — y - w 

(T3) = W++ - w+- - W-+ +w — , 



(2) 



such that etc. can be found from the average values of 
the three observables {T„i). An eigenstate of the observable 
A has a fixed value for all measurements, or a vanishing 
dispersion (A^) - (A)^ = , (A) = ±1. For (Ti) =1 the 
probability pr for all classical states for which Ar — —1 
must vanish. Such a state leads to w ^ = w = 0. 

Pure and mixed states 

Let us now specify our subsystem. For the manifold 
of all pure states of the subsystem we choose the ho- 
mogeneous space SU{A)/ SU{3) x U{1). We parameter- 
ize the embedding space by the 15 components fk 
of a vector {k — 1 . . . 15). It is normalized according to 
J2k fk — ^'^^ obeys eight additional constraints that 
reduce the independent coordinates to six, as appropri- 
ate for the dimension of the complex projective space 
= SU{4)/SU{^) X U{1). (For a geometrical discussion 
of the complex projective spaces for pure states in quan- 
tum mechanics cf. ref. [ll|.) An easy way to obtain the 
constraints for fk employs a hermitean 4x4 matrix p. 



fkLk) , fk = tT{pLk 



(3) 



(Summation over repeated indices is always implied.) Here 
Lk are fifteen hermitean 4x4 matrices obeying 

Ll = l , trLfc - , tiiLkLi) = iSki. (4) 

They read explicitly (with the Pauli 2x2 matrices) 



L3 



diag(l, 1, -1,-1) , L2 - diag(l, -1, 1,-1) , 

n, \ 
0, Tl ) ' 



diag(l, -1,-1,1) , L4 



(5) 



T2, 
0, T2 



n, \ , f T2, 



with Lg, Lg, Z/10, ill obtained from (L^, L^, Lq, Ly) by ex- 
changing the second and third rows and columns, and 
Li2, Li3, Li4, _Li5 similarly by exchange of the second and 
fourth rows and columns. The matrix p parameterizes the 
homogeneous space SU{4:)/ SU(3) x U{1) if it obeys 

p = , [/[/t ^u'fU , pi^ diag(l, 0, 0, 0), (6) 

for some appropriate unitary matrix U . This implies 



y2 Pac 
^ — ^ a. 



trp = 1. 



(7) 



The observables ^1,2,3 are specified by {Tm){fk) = 
fm , 771 = 1, 2, 3, which is equivalent to the specification of 
w'^"^\fk) in eq. ([1]). Already at this stage we get a glance 
on the possibility of entanglement, since states with /i = 
/2 = , /3 = -1 will lead to (Ti) = (T2) = 0, (T3) = -1, 
and therefore to a correlation for opposite values of bit 1 

and bit 2, w++ = w = , = w |- = i. We label 

these two-level observables by a real vector with compo- 



nents Cfc, with ek{T„ 



1...3 , k 



1. 



15. 



The mean value of Tm in a given state fk can then be 
written in the form (with e^™"* = ek{Tm)) 



{Tm){fk) ^ fke'r\ 



(8) 



We next represent an observable A{ek), labeled by e^, in 
terms of a hermitean operator 



A = ekLk , ek{A) ^ itr(iLfc), 



(9) 



with J2k^t ~ ^ ^'^'^ A^ = 1. In this language one finds 
i'-^m) — tr(T'^p), = 

The states labeled by fk will correspond to pure states 
in the four-state quantum system. Mixed quantum states 
obtain, for example, by the replacement fk — > Pk with 



Pk 



1/2 



fk, 



and purity P obeying 



P = PkPk < 3. 



(10) 



(11) 



The state of the subsystem is now determined by the fifteen 
numbers pk and we assume for m = 1, 2, 3 that (Tm) = pm- 
For given pk the state of the subsystem can be described 
by a density matrix p, 



P= ^(1 + Pfeife)- 



(12) 



In terms of p one finds the familiar quantum law for expec- 
tation values 



{A) = tr(ip) 



(13) 



for every observable whose expectation value obeys {A) = 
PkSk, as for the "one bit observables" T,n- 

We also can use p in order to formulate the general con- 
dition that the pk have to obey if they describe a four-state 
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quantum system. An arbitrary hermitean p can be diago- 
nalized by a suitable unitary transformation 



U diag (pi 



)C/t 



and we require that p is a positive matrix, 
0<pa<l. 



(14) 



(15) 



The relations (fT4|) . (fT5|) . together with the definition 
(fT2|l . may be called "positivity constraint" or "purity con- 
straint" . A subsystem characterized by expectation val- 
ues pk can behave as a quantum system only if the fifteen 
numbers pk obey the purity constraint. In turn, the purity 
constraint ensures that the diagonal elements of the den- 
sity matrix are positive and bounded by one, < paa < 1- 
For pure states one has P = 3 and p coincides with p ([7]). 
More generally, one finds 



P ^ PkPk = 4trp2 - 1 



1 < 3. 



(16) 



The purity is a statistical quantity characterizing the fluc- 
tuations in the subsystem. We will see that it limits the 
number of independent system observables which can have 
sharp values simultaneously, motivating the name "purity 
constraint" . 

The description of pure states in terms of wave functions 
ij^a , V'^V' = 1; can be obtained from the density matrix in 

a standard way, Pap = Ipa^p , Ipa = Uaf3{i'l)f3 , (V'm)a = 

6ma- This expresses the fk as a quadratic form in the 
complex four- vector ipa, 



(17) 



and shows directly that only six components of fk are in- 
dependent. The quantum mechanical wave function i/; ap- 
pears here as a convenient way to parameterize the man- 
ifold of pure states. The expectation value (A) = fkSk 
translates into the quantum mechanical law 



(A) = V'^iV'- 



(18) 



Assume now that the time evolution of the classical prob- 
ability distribution {pr} changes a pure state of the subsys- 
tem into another pure state. It therefore corresponds to a 
rotation in the space of the (constrained) vectors fk which 
parameterize S';7(4)/S't/(3) X C/(l). According to eq. ^ 
this is represented in the quantum language by a unitary 
evolution of the pure state density matrix 



pit) = U{t,t')pit')U\t,t'). 



(19) 



The Schrodinger equation for the wave function can be eas- 
ily derived from eq. The unitary time evolution 
is a characteristic feature for the isolation of the subsystem. 
We may assume that it also holds for mixed states, i.e. all 
density matrices obeying the purity constraint. The con- 
sequences of this setting are simple but striking: the time 
evolution of expectation values of the system observables 
precisely follows the time evolution in quantum mechanics 
for an appropriate Hamiltonian H = i{dtU {t^t' fjU'^ {t,t'). 



If H is independent of p^, the quantum time evolution is 
linear. 

Entangled states 

Let us next look at a typical entangled quantum state. 
We consider the wave functions 



1 

71 



{^2 ± -03 



(20) 



with associated pure state density matrices p±. These 
states are eigenstates to T3 with eigenvalue —1. Writing 



P± = J(l-i3±(il2-il4)), 



(21) 



we infer /s = -1 , /12 = ±1 , /14 = Tl, and aU other fk 
vanishing. Thus (Ti) = (T2) = implies that the values of 
bit 1 and bit 2 are randomly distributed, with equal prob- 
abilities to find -1-1 or —1. Nevertheless due to (T3) = — 1, 
the product of both bits has a fixed value. Whenever bit 
1 is measured to be positive, one is certain that a mea- 
surement of bit 2 yields a negative value, and vice versa. 
The two bits are maximally anticorrelated. The situation 
is the same for "spin observables in other directions" that 
we will discuss later in this note. Obviously, an entangled 
state can be realized by a subsystem with the appropriate 
values Pk = fk, and its time evolution will be unitary if the 
purity is conserved. 

At this stage we may give an explicit example for a clas- 
sical statistical ensemble that realizes such an entangled 
state. Let us label the classical states as r = ({CTfc},C), 
with {cTfc} an ordered sequence of fifteen discrete variables 
(Tfc = ±1 which accounts for 2^^ different possibilities. The 
number of classical states is typically higher than 2^^ - this 
is taken into account by the additional index C,. We identify 
the expectation values which specify the state of subsystem 
as 



Pk 



^ akp{{(Tk},C), 

({'^fc}:C) 



(22) 



with Pk = fk for pure states. In particular, the two bits or 
two spin variables Ti , T2 are associated with the classical 
observables cti and (T2, while specifies the measurement 
correlation of pairs of measurements for Ti and i.e. 

{T2Tl)m = P3 = (cts)- 

The probability distributions {pr} which describe an iso- 
lated subsystem (these are not the most general allowed 
probability distributions) can be written as 



p{Wk},C) 

PsiWk}) 
^ (JkSpe{{(Jk},C) 



Ps{Wk})Ps{0 +Spe{{ak}X), 

2-15[|(l+CTfcPfc), 

fe 

, 5]p,(C) = 1. (23) 

C 



The information contained in the part of the probability 
distribution Spe describes properties of the environment, 
but is not needed for the computation of the state of the 
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subsystem. We observe that the classical correlation func- 
tions as {(JkCTi) typically depend on Spe- They cannot be 
computed from the information available for the subsys- 
tem, i.e. in terms of pk- The subsystem is described by 
"incomplete statistics" in this sense. We will argue be- 
low that sequences of measurements in the subsystem are 
described by a different conditional correlation function. 

Finally, the unitary time evolution characteristic for the 
isolated subsystem is realized by an evolution oips{pk) ac- 
cording to the unitary evolution of the pk corresponding to 
eq. p^ . while the evolution of Spf. is arbitrary as long as 
the constraint is respected. The classical probability 
distribution for the entangled state ([2T|) obtains by insert- 
ing va.ps{pk) the appropriate values pk = fk- Since an arbi- 
trary unitary evolution can be described by an appropriate 
evolution of the classical probability distribution {pt}, we 
can realize "quantum operations" as the "CNOT-gate" for 
the two bits by a suitable time evolution of a classical sta- 
tistical ensemble. It is sufficient that after a given time 
interval the unitary matrix takes the form 



U 



( 1,0,0,0' 
0,1,0,0 
0,0,0,1 

V 0,0,1,0, 



(24) 



Interference 

Other interesting quantum phenomena are the superpo- 
sition of states and interference. Just as entanglement, 
they can be described within our classical statistical en- 
semble. We concentrate on pure states p^ — fk- Consider 
two pure quantum states that evolve in time according to 

i^a = -i=(^i-f ^2)e-^""* , Vfc = -^(V'i-V'2)e-'^"''*. (25) 

The corresponding density matrices are time independent. 



Pa,b^il + L,±Li±Le)/4. 



(26) 



Both states describe an eigenstate of the first bit, (Ti) = 1, 
whereas the second bit is randomly distributed, (T2) = 0. 
Due to the time dependent phase, the interference can be 
positive or negative for the superposition of both states, 
tp = ■^{''Pa +'^Pb)- For such a superposed state a quantum 
mechanical computation leads to a characteristic oscillation 
of {T2), 

{T2) = ^^L2^ = COs(Ai) , A = LJa - C^b, (27) 

as known from the oscillation of a spin in the z-direction for 
a superposition of spin-eigenstates in the x-direction, with 
different energies for the positive and negative Sx eigen- 
values. A rotation of the pk = fk, describing the time 
evolution of the classical probability distribution {pr} ac- 
cording to eq. with f2 — fs — cos(Ai), reproduces 
the "interference pattern" ((27|) . An evolution law dtf2 = 

A/5 , dtf5 = -A/2 , /3 = /2, /y = /s , /i = 1 , fk = 

otherwise, has indeed solutions leading to a density matrix 
which corresponds to the superposed state '0, 



P 



i 1 1 +Li -Fcos( At) {L2+L3) - sin( At) (L5 -HL7) | . 



(28) 



Bosons and fermions 

Furthermore, our classical statistical ensemble can de- 
scribe identical bosons or fermions. We may identify the 
two bits with two particles that can have spin up or down, 

^1 = I TT) , ^2 = I U) , ^3 = lit) , V'4 = I U). (29) 

If the particles are identical, no distinction between bit 1 
and bit 2 should be possible. This requires that the sys- 
tem must be symmetric under the exchange of the two bits, 
imposing restrictions on the allowed probability distribu- 
tions {pt}- The symmetry transformation corresponds to 
an exchange of the second and third rows and columns 
of p. On the level of the fk this amounts to a mapping 
fk^ f'k ■■ /i ^ /2 , /4 ^ /s , /s ^ /g , /e ^ /lo , /? ^ 
/ii , /i3 ^ /i5, while /3,/i2 and /u remain invariant. 
Allowed probability distributions ([^5]) must obey ps(/fc) = 
Psifk)- In particular, the allowed pure states are restricted 
by A = /2 , /4 = /s , /s = /g , /e = fio , /t = fii, and 

/l3 = /l5- 

Consider the pure states ?/'+ and ip- in eq. (^0)1 . For 
both states the density matrix p± is compatible with the 
symmetry. This does not hold for the density matrices 
corresponding to the states -02 or ip^. In fact, linear super- 
positions of ■0+ and ijj- are forbidden by the symmetry, a 
pure state aip- + bijj^ must have a = or 6 = 0. The sym- 
metry requirement acts as a "superselection rule" for the 
allowed pure states or density matrices. For an arbitrary 
state vector aip- + bipj^ + cijji + (iV'4 symmetry of p re- 
quires either a = 0or6 = c = rf = 0. We observe that ip- 
switches sign under the symmetry transformation as char- 
acteristic for a state consisting of two identical fermions. 
In contrast, the boson wave function ■0 = bip^ + cipi + dip2 
is invariant under the "particle exchange symmetry" . The 
time evolution for "identical bits" must be consistent with 
the exchange symmetry. It corresponds to SU{?>) rotations 
in the space of bosons, while the fermion state has a time 
independent density matrix. 

System observables 

So far wc have only used the mutually commuting opera- 
tors Li, L2, L3. One may ask if other operators of quantum 
mechanics, as arbitrary L^, can correspond to observables 
in our classical statistical setting. It is straightforward to 
find probabilistic two level observables Tk corresponding 
to arbitrary Lk- We require (Tk) = pk and associate Tk 
with the classical observables ak in eq. (|22p . Associating 
to every quantum operator the vector according to eq. 
©, we specify (Tk) — Pk^k, leading directly to the quan- 
tum rule Uni) for {Tk). We can obviously find eigenstates 
for arbitrary Tk, as pa in eq. (j26p . which is an eigenstate 
to T4 with eigenvalue -1-1 and similar for Tg. In fact, the 
operators Li, L4, Lg mutually commute such that there ex- 
ist simultaneous eigenstates for all three observables. On 
the other hand, T2 cannot have a sharp value in this state 
- L2 does not commute with L4 and Lg. The quantum 
mechanical uncertainty principle for non-commuting op- 
erators is directly implemented in our classical statistical 
setting. If we associate (Li, L2) to spins in the z-direction, 
it is straightforward to associate (^8,^4) and {Lq,Lz) to 
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spins in the x- and y-directions, respectively. We may- 
consider the fifteen observables Tj. as basis observables. A 
sharp value for a basis observable requires pk = ±1, and 
we see how the purity constraint (|16p limits the number of 
simultaneously sharp basis observables to at most three. 

While the observables Tk are again two-level observables 
with spectrum -1-1 and —1 (corresponding to = 1), the 
possible classical observables associated to arbitrary her- 
mitean operators can have a spectrum of possible values 
for individual measurements with up to four different real 
values. Furthermore, we may also wish to consider rotated 
two level observables in arbitrary directions. The realiza- 
tion of such a larger set of observables of the subsystem 
as classical observables with fixed values in every state r 
of the classical ensemble requires infinitely many classical 
states r Q (using the index ( in eq. (|22p). One finds 
that arbitrary observables of the four-state quantum sys- 
tem can find a classical realization in this way. To all such 
system observables one can associate a vector such that 
{A{ek)) — PkSk- Therefore the "quantum rule" for expec- 
tation values, (A(efe)) = tr(A(efc)p), with A{ek) = CkLk, 
holds for all system observables. In particular, the prob- 
abilities to find for arbitrary two-level observables A{ek) 
(represented by operators obeying A^{ek) — 1) the values 
+1 or -1 obey w±{pk; e^) = i(l ± PfcCfc). 

One may wonder why the "system observables", whose 
expectation values obey the simple rule (A) = pkCk, play 
such an important role for the description of subsystems 
which show a quantum character. We invoke here the 
"principle of equivalent state and observable transforma- 
tions" (PESOT). This states that it should be possible to 
undo any time evolution of the state by a corresponding 
time evolution of the observables. In the quantum me- 
chanical language, this amounts to the well known equiv- 
alence between the Schrodinger and Heisenberg pictures - 
the time evolution can be either described by the evolution 
of the state or the evolution of the operators. According 
to PESOT one can always define an appropriate observ- 
able transformation such that (A) remains invariant under 
a change of . Since the pk transform under SU (4) as 
the adjoint representation (with additional nonlinear con- 
straints (Uni), (fTTjl ). one needs another adjoint repre- 
sentation characterizing the transformation of the observ- 
ables, namely e^, in order to construct the invariant scalar 
product. Together with linearity in (or pk) this fixes 
(A) = CkPk- 

Conditional correlations 

Beyond a rule for the computation of expectation values 
a probabilistic description of Nature has to specify a rule 
how the probability for sequences of events or sequences 
of measurements is computed. For a sequence of measure- 
ments of first B and then A (which may be associated to a 
sequence of two events) we need the conditional probability 
to find one of the possible measurement values for A under 
the condition that a certain value of B has been measured. 
In general, conditional probabilities depend on the precise 
way how measurements are performed, since the measure- 
ment of B changes the subsystem and this change has to 
be specified. For a very large number of effective degrees 



of freedom, as relevant for classical thermodynamics, the 
influence of the measurement of B on the system is often 
neglected. This is not valid for our description of subsys- 
tems. 

We denote the conditional probability to find a value 
e for bit 2 if bit 1 has been measured to have a value 
7 by p{e;j). For the entangled state ((20|) . (|2ip it obeys 
p(l; 1) = p{-l; -1) = , p(l; -1) = p{-l; 1) = 1. In our 
classical statistical setting we describe the product of two 
measurements of observables A and i? by a "measurement 
correlation" {AB)m.- For the two bits Ti and T2 it is given 
by ((Ti) = wi+ - 

(72 Ti)™ = p(l; l)wi+ +p(-l; 

-p{l--l)wi^-p{-l;l)wi+, (30) 

and we can identify the probabilities in eq. ^ as w++ — 

p(l; l)u)i+ , — p(l; —l)wi- etc.. In turn, for wi± 7^ 

the conditional probabilities p(e; 7) can be computed from 
{T2Ti),n = Pz and (Ti) = pi , {T2) = P2, using p(l;7) + 
p(— 1;7) — 1. Here we assume that for the special case 
of two "comeasurable" bits the expectation value {T2) is 
not influenced by the measurement of Ti, such that {T2) — 
p{l\ l)wi+ + p(l; -l)wi- - p(-l; l)wi+ - p(-l; -l)wi-. 
This will not hold for general pairs of observables, in par- 
ticular not for pairs which will be represented by non- 
commuting quantum operators. 

In general, we require that the outcome of a sequence of 
measurements which respect the isolation of the subsystem 
should be predictable for a given state of the subsystem, 
and not involve details of the environment. In other words, 
the measurement correlation {AB)„i for a pair of two "sys- 
tem observables" should be computable in terms of pk and 
be independent of Sp^. Furthermore, we require that the 
repetition of a measurement yields an identical result, and 
that good measurements keep the maximum amount of in- 
formation about the subsystem which is compatible with 
the first two requirements. One can then show Q that the 
measurement correlation equals the quantum correlation, 
as expressed by the expectation values for the anticommu- 
tators of operators, 

(AB)™-itr({i,B}p). (31) 

The probabilities for the outcome of two measurements fol- 
low therefore exactly the quantum mechanical laws and can 
violate Bell's inequalities [T^ . 

In contrast, Bells inequalities apply to classical correla- 
tion functions for classical observables [3]. However, as we 
have seen before, the classical correlations depend on Sp^. 
They cannot be used for measurements which maintain the 
isolation of the subsystem and only measure properties of 
the subsystem. The classical correlations are appropriate 
only for measurements of properties of the subsystem to- 
gether with properties of the environment, in which we are 
not interested. Concerning locality and realism, our clas- 
sical statistical setting shares the same properties as other 
classical statistical systems. Correlations can be non-local 
- it is sufficient that they have been prepared in the past 
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by causal processes. Correlated systems have to be con- 
sidered as a whole - parts of such a system cannot be as- 
sociated to independent entities, even if signals cannot be 
exchanged any longer between the parts. Correlations are 
part of reality on the same footing as expectation values 
of observables. The reason for the violation of Bell's in- 
equalities by the measurement correlations can be rooted 
in the incompleteness of the statistical description of the 
subsystem even for ensembles where the description of the 
subsystem together with its environment is complete. 

In conclusion, we have derived four-state quantum me- 
chanics from classical statistical mechanics. Every quan- 
tum system with arbitrary Hamiltonian and observables 
finds an equivalent description in terms of a classical sta- 
tistical ensemble with an appropriate time evolution of the 
probability distribution. In particular, we have presented 
an explicit description of entanglement and interference, as 
well as fermions and bosons, within a classical ensemble. 

Our approach is not a simple rewriting of the quantum 
mechanical laws. Quantum statistics appears as a particu- 



lar case within a much larger class of probabilistic theories 
for subsystems. The restriction to a manifold of pure states 
SU{M)/SU{M - 1) X U{\), the unitary time evolution of 
mixed states, or the principle of equivalent state and ob- 
servable transformations could be abandoned, leading to 
statistical models different from quantum mechanics. For 
example, a decrease of the purity describes decoherence 
jl4] , while in the opposite case an increase of purity results 
in the approach towards a pure state or "syncoherence" 
Q. Non-linear versions of quantum mechanics could be 
described if the unitary matrix U{t, t') in eq. (fT9|) depends 
on pfe. Our approach therefore provides a framework to 
test quantum mechanics against "neighboring" statistical 
ensembles where some of the properties leading to quantum 
mechanics are not obeyed exactly. Besides the conceptual 
and perhaps philosophical impact of a classical statistical 
description of quantum mechanics, it remains to be seen 
if our findings can also be of practical use, for example 
by classical computations of quantum processes or steps in 
quantum computing. 
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